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I.   INTRODUCTION 

Numerous  methods  of  estimating  tropospheric  propagation 
path  loss  are  currently  found  in  the  literature.   These 
include  both  graphical  and  computer  techniques,  of  which  the 
graphical  methods  are  most  widely  presented.   In  conducting 
this  particular  study  a  program  dealing  with  tropospheric 
path  loss  prediction  (TR0P0PL0T)  was  chosen  as  a  computer 
method  of  analysis.   This  program  was  originally  published  by 
ESSA  in  a  technical  report  [Ref.  6]  and  was  subsequently 
modified  for  use  at  the  Naval  Postgraduate  School  [Ref.  7]. 
Tests  were  conducted  by  Longley  and  Reasoner  [Ref.  8],  in 
which  the  computer-predicted  results  of  the  program  were 
compared  with  empirical  data  taken  over  a  number  of  different 
paths.   The  results  showed  that  TR0P0PL0T  can  provide 
reasonably  accurate  predictions  of  path  loss  within  certain 
constraints. 

The  paths  selected  for  analysis  in  this  study  are  located 
in  the  Puget  Sound,  Washington  area.   They  comprise  both 
existing  and  proposed  communications  links  connecting  the 
underwater  range  facilities  operated  by  the  Naval  Torpedo 
Station  at  Keyport,  Washington.   The  range  operation  centers 
are  currently  located  at  Winchelsea  Island,  British  Columbia 
and  on  Zelatched  Point  near  Dabob  Bay.   The  present  communica- 
tions requirements  include  on-range  communications  with  both 
range  control  vessels  and  the  submarines  and  aircraft  conduc- 
ting tests,  the  capability  to  monitor  telemetry  information, 


12 


and  the  maintenance  of  telephone  links  with  other  sites. 
Proposals  for  expansion  of  the  range  facilities  as  the 
TRIDENT  program  progresses  include  possible  at-sea  range 
facilities  and  digital  data  transmission  to  a  central  computer 
processing  center  located  at  Bangor  or  Keyport.   These  pro- 
posals would  increase  the  bandwidth  requirements  of  the  links 
and  thus  necessitate  a  change  from  the  present  VHF  frequency 
range  to  an  S-band  (2  GHz)  line  of  sight  link. 

The  terminal  and  repeater  sites  are  shown  in  Figure  1. 
The  existing  VHF  link  consists  of  paths  from  Winchelsea  to 
Lookout  Mt.  (repeater)  to  Bangor  and  from  Zelatched  Point  to 
Bangor.   The  proposed  link  at  2  GHz  runs  from  Makah  to  Striped 
Peak  (repeater)  to  Mt.  Constitution  (repeater)  to  Gold  Mt. 
(repeater)  to  Bangor  or  Keyport,  as  well  as  from  Zelatched 
Point  to  Bangor  or  Keyport.   Each  of  these  paths  is  computer- 
analyzed  in  Section  II. 

In  addition  to  computing  path  loss,  the  historical 
meteorological  data  for  this  area  was  examined  in  an  attempt 
to  characterize  the  "typical"  propagation  conditions  found  in 
this  vicinity.   This  was  accomplished  by   means  of  the 
Integrated  Refractive  Effects  Prediction  System  (IREPS)  which 
provided  output  based  on  data  taken  from  upper  air  sounding 
stations  in  the  Puget  Sound  area.   This  information  could  also 
be  useful  in  the  determination  of  link  service  probability. 
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Figure  1 
NAVTORPSTA  Keyport  Ranges  and  Vicinity 
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A  listing  of  a  computer  program  designed  to  plot  path 
profiles  and  to  aid  in  the  determination  of  several  input 
variables  is  included  in  the  appendices,  as  is  an  explanation 
of  the  modifications  required  in  order  to  use  detailed  terrain 
profile  information  in  the  TR0P0PL0T  program. 
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II.   COMPUTER-AIDED  PATH  LOSS  PREDICTION, 

A.    PROGRAM  DESCRIPTION 

TROPOPLOT  provides  as  an  output  a  measure  of  the  long- 
term  attenuation  over  a  given  path.   The  program  was  primarily 
designed  for  use  in  cases  where  detailed  terrain  information 
was  not  available  but  can  be  used  accurately  with  user- 
supplied  data  after  a  slight  modification.   The  required 
input  parameters  include  frequency,  antenna  heights,  antenna 
gains,  line  losses,  receiver  sensitivity,  path  length, 
polarization,  transmitter  power  out,  surface  ref racti vi ty  , 
conductivity,  permittivity,  and  a  measure  of  the  terrain 
roughness  Ah.   These  parameters  are  for  the  most  part  well 
defined  in  Refs.  6  and  7  and  are  easily  determined  except  in 
the  case  of  user-supplied  terrain  profile  information.   For 
this  sort  of  input  there  is  some  ambiguity  concerning  the 
selection  of  values  for  a,  e,  Ah,  and  the  effective  antenna 
heights. 

In  considering  a  path  which  includes  both  poor  ground 
surfaces  (e  =  4,  a=.001)  and  sea  water  surface  (e=81,  a=5)  the 
particular  value  to  use  as  input  is  not  specified,  although 
a  completely  oversea  path  requires  an  adjustment  to  the 
program  not  available  in  the  version  used  for  this  study. 
Thus  the  values  of  e  and  a    chosen  were  based  on  the  corres- 
ponding value  for  that  portion  of  the  terrain  that  constituted 
the  largest  portion  of  the  dominant  reflecting  plane  between 
the  transmitter  and  receiver. 
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The  path  geometry  considered  is  shown  in  Figure  2  and 

indicates  h  -,  9  as  the  structural  transmitter  (receiver) 
g  I  ,  c 

height  above  ground.   This  input  parameter  is  converted  into 

an  effective  height,  h  ,  «  b^  tne  program.   It  should  be 

noted  that  Longley  and  Rice  consider  two  cases:  one  for 

random  antenna  siting  in  which  structural  heights,  h  ,  0  are 

gl  ,2 

considered  equal  to  the  effective  heights,  h  ,  ? ;  and  another 
in  which  the  sites  are  carefully  selected,  as  in  the  case  of 
radio  relay  links,  and  the  effective  heights  are  larger  than 
the  structural  heights  as  shown  below: 

h  i  n    ~    h,  0  +  kexp(-2h  ,  0/Ah)  meters 
el, 2    g 1  ,2      rx   g 1  ,2 

The  variable  k  is  considered  to  have  a  maximum  value  of 
50  as  determined  by  the  author's  study  of  varied  terrain 
conditions.   The  method  presented  below  for  the  computation 
of  k  is  valid  only  for  antenna  heights  less  than  or  equal  to 
10  meters.   In  Callaghan's  version  [Ref.  7]  this  determination 
was  : 


k  =  1  +  4  s  i  n  (  tt  h  .  o/10)  for  0  <  h  1,2  <  5 

gl  ,2'         -  g    - 


=  5 


otherwi  se 


Since  most  of  the  antennas  in  the  links  under  study  had 

structural  heights  in  the  neighborhood  of  50  meters,  it  is 

obvious  that  predictions  based  on  strictly  h  1  7   would  in  this 

g  i  ,  c 

case,  produce  an  erroneous  (higher  than  normal)  prediction  of 
path  loss.   This  factor  is  recognized  by  the  authors  as  an 
area  requiring  better  definition  and  as  a  primary  source  of 
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Figure   2 
Path  Geometry 
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prediction  error.   For  this  reason  the  modifications  to 
TR0P0PL0T  which  are  detailed  in  Appendix  B  were  necessary. 
The  effect  of  varied  antenna  heights  in  the  input  of  the 
unmodified  program  is  shown  in  the  output  data  presented 
later  in  this  section. 

The  description  of  the  terrain  is  accomplished  by  a 
statistical  quantity,  Ah,  which  is  defined  in  Ref.  6  as  "the 
asymptotic  value  of  the  interdecile  range,  Ah(d),  of  terrain 
heights  above  and  below  a  straight  line  fitted  to  elevations 
above  sea  level.   The  parameter  Ah(d)  is  calculated  at  fixed 
distances  and  its  median  value  usually  increases  with  path 
length  to  Ah."   For  a  single  path  profile  these  definitions 
are  not  adequate.   The  asymptotic  value,  Ah,  used  in  this 
study  as  determined  by  taking  the  interdecile  range  of  the 
difference  between  a  straight  line  fitted  to  the  path  profile 
points  and  the  path  profile  points  themselves.   The  results 
of  this  estimation  of  Ah  correlate  favorably  with  those 
values  listed  in  Table  1  [Ref.  7]  which  shows  estimated  values 
of  Ah  for  particular  types  of  terrain.   Figure  3  shows  an 
example  of  the  straight  line  fit  to  the  terrain  profile  between 
Lookout  Mt.  and  Bangor.   The  details  of  the  computation  are 
found  in  Appendix  B. 

Having  thus  considered  the  variations  in  the  input 
parameters  the  program  constraints  should  be  noted.   TR0P0PL0T 
is  designed  for  use  within  the  following  constraints: 
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Parameter 

frequency  (f) 

antenna  height  ( h  ,  0) 

g  1  ,2' 

distance  (dist) 

surface  refractivity  (N  ) 


Range 
20  to  40,000  MHz 
0.5  to  3000  m 
1  to  2000  km 
250  to  400  N-units 


1  ,2 


The  antenna  siting  is  subject  to  the  following  conditions 

1.  The  angle  of  elevation,  e  1    n ,    of  each  horizon 

3  el  ,2 

ray  from  the  horizontal  should  not  exceed  12°. 

2.  The  distance  from  each  antenna  to  its  horizon 
(d.,  2)  should  not  be  less  than  1/10  or  more 
than  3  times  the  corresponding  smooth  earth 
distance  (d, s ) • 

These  limitations  can  be  ignored  if  the  values  of  d. 

and  e  ,  «  are  entered  directly  into  the  program.   As  is  the 
el  ,2  ■*  r   3 

case  with  Ah  and  h  ,  0  the  method  of  accomplishing  this  is 

gl  ,2 

contained  in  Appendix  B. 

Having  resolved,  or  at  least  indentified,  the  potential 
ambiguities  in  the  input  parameters,  it  is  useful  to  briefly 
describe  the  output  format  of  the  program.   As  shown  in 
Figs.  4-7,  the  output  of  TR0P0PL0T  consists  of  two  tables 
containing  transmission  loss  versus  distance  and  signal 
strength  versus  distance,  a  plot  of  the  transmission  loss 
information  and  a  printout  of  the  calculated  and  input  values 
of  many  of  the  program  parameters.   Some  of  the  quantities  in 
this  figure  may  require  some  explanation  as  to  their  meaning 
or  derivation.   These  include  the  variables  AE,  AES,  AED,  MS, 
~*  MD,  ADX,  DX,  Kl  ,  K2  ,  and  ALS. 
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In  the  process  of  describing  these  parameters  a  brief 

outline  of  how  TR0P0PL0T  determines  path  loss  is  presented. 

A  more  detailed  exposition  can  be  found  in  Ref.  6.   The 

median  reference  value  of  attenuation  below  free  space,  A   , 

r      cr 

is  computed  first.   The  reference  value  of  transmission  loss, 

L   ,  then  becomes  the  sum  of  the  free  space  attenuation,  L,f, 

and  the  reference  attenuation,  A   , 

cr 


L„   =  L.  ,    +    A    dB. 
cr    bf    cr 


where  the  free  space  loss  is  defined  as, 

Lbf  =  32.45  +  20  log1Q  (f  in  MHz)  +  20  log1Q(d  in  km) 

The  reference  attenuation,  A   ,  is  determined  using  one 
of  three  subroutines  depending  on  the  particular  mechanism 
of  propagation;  LOS  for  line  of  sight  modes,  DIFF  for 
diffraction,  and  SCATT  for  tropo-scatter .   Two  ray  optics  is 
used  to  compute  line  of  sight  paths  while  diffraction  paths 
are  assumed  to  be  over  a  double  knife  edge  and  a  value  of 
diffraction  attenuation  below  free  space,  A,,  is  computed. 
In  the  case  of  scattering  the  attenuation  variable  is  designa- 
ted as  A  .   The  reference  attenuation  A   is  determined  by 
s  cr  J 

the  smaller  value  of  A,  and  A  .   In  the  line  of  sight  case 

d      s  3 

the  attenuation  is  calculated  for  2  values  of  distance 
(d  ,  d ,  )  for  which  line  of  sight  propagation  is  valid  to 
produce  corresponding  values  of  attenuation  A   and  A-,.   The 
diffraction  attenuation,  A.,  is  computed  at  the  distance 
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d.   and  together  with  A   and  A-,  is  used  to  determine  the 
slopes,  k,  and  k?,  of  a  smooth  curve  of  the  reference  attenua 
tion  versus  distance  over  the  range  1  <_  d  <_  d 

Acr  =  Ao  +  kl(d_do)  +  k2  log10(d/do)  dB' 


The  output  parameter  A   is  defined  as 

e 


A   =  A   -  k,d   -  k0  log,nd    (i.e.,  d  =  l  ) 
e    o    lo    2   31 0  o   x         ; 


so  for  1  <  d  <  d 


ls 


A   =  A  +  k, d  +  k0  log,nd 
cr    e    1     2   31 0 

When  the  diffraction  attenuation  is  computed  the  output 

variables  A   and  m,  come  into  consideration.   The  diffraction 
ed      d 

attenuation  is  computed  as  the  weighted  average  of  the 
diffraction  attenuation  over  smooth  earth,  A  ,  and  the  attenu- 
ation over  a  double  knife  edge  surface,  A,  ,  where 


A  ,  =  (1 -w)A.  +  wA   dB 
d       J    k     r 


A  description  of  the  method  of  weighting  these  estimators- 

and  calculating  A,  and  A.  is  contained  in  Annex  3  to  Ref.  6. 
3   d      k 

The  diffraction  attenuation,  A,,  is  determined  at  two 

distances,  d   and  d„ ,  in  the  far  diffraction  region,  and  a 
3      4 

straight  line  through  the  points  (A-,d»)  and  ( A  -  ,  d .  )  is 
defined  as 


md  =  (A4  "  A3)/(d4  "  V  dB/km    (slope) 
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and 


ed    to    4    d  4 


( i  intercept) 


where  A,   is  a  "clutter  factor"  (<  15  dB). 
to  — 

If  the  scatter  attenuation  is  less  than  the  diffraction 

attenuation,  as  is  sometimes  the  case  in  trans -hori zon  paths 

where  the  distanced  or  the  angular  distance  9  is  large,  then 

A   =  A  .   When  the  product  of  d  in  kilometers  and  9  in 
cr    s  r 

radians  is  greater  than  0.5,  A   is  computed  at  large  distances, 
d5  and  dfi,  and  a  straight  line  through  the  points  ( A  j-  ,  d  j- )  and 
(Ag,dg)  is  defined  as  follows: 

ms  =  <A6  "  V<d6  "  V 


and 


Aes  =  A5  "  m5d5 


As  before  the  reference  attenuation  then  becomes 


A   =  A   =  A   +  m  d 
cr    s    es     s 


d  >  d 
—  x 


The  quantity  d   is  defined  as  the  distance  where  the  scatter 
attenuation  is  equal  to  the  diffraction  attenuation. 

It  can  be  seen  from  the  above  description  and  Ref.  6  that 
the  actual  determination  of  the  output  quantities  is  somewhat 
complex  and  for  most  applications  the  information  contained 
in  the  graphs  and  tables  describing  path  loss  and  signal 
strength  is  sufficient.   Consequently  the  data  presented  in 
the  next  sub-section  will  primarily  consist  of  derivations 
from  those  portions  of  the  output.   If  the  quantities  in 
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Fig.  6  are  important  to  the  user  it  should  be  noted  that  for 
certain  input  combinations  the  format  statements  contained 
in  the  TR0P0PL0T  program  will  produce  an  all  asterisk, 
11*******1'^  printout  indicating  that  the  format  statement 
governing  that  particular  variable  requires  modification. 

B.    ACTUAL  PATH  LOSS  MEASUREMENTS 
1 .    Procedure 

In  order  to  obtain  the  proper  input  parameters  for 
entry  into  TR0P0PL0T  certain  preliminary  measurements  and 
calculations  were  required.   The  first  step  was  to  obtain 
charts  of  the  path  area   and  from  these  determine  the  terrain 
profile  over  a  great  circle  path  between  transmitter  and 
receiver.   The  great  circle  path  was  approximated  by  a  rhumb 
line  for  path  lengths  less  than  70  kilometers  and  the  avail- 
ability of  7V  topographic  charts  for  the  longer  paths  enabled 
the  use  of  straight  lines  over  the  chart  area.   This  profile 
information  was  plotted  on  a  curved  earth's  surface  using 
the  plotting  program  described  in  Appendix  A.   The  value  of 
the  surface  refractivity  used  in  determining  the  effective 
earth's  radius  was  the  same  as  that  used  in  the  input  to  the 
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To  obtain  a  catalogue  for  U.S.  Charts  write  to 

Denver  Distribution  Section 

U.S.  Geological  Survey 
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program  and  can  readily  be  obtained  from  Ref.  9  or  from 

actual  meteorological  data.   If  the  charts  of  Bean  e_t  al  , 

(1960)  are  used  this  value  must  be  converted  from  N   to  N   by 

os 


the  equation 


N$  =  NQ  exp(-0.1057  h  ) 


The  value  of  h   used  depends  on  several  factors  most  important 
of  which  is  the  mode  of  propagation.   For  most  of  the  paths 

considered  herein  the  value  of  h   chosen  was  determined  by  the 

s  J 

elevation  of  the  lowest  antenna.   For  trans-horizon  paths  a 
mean  N   is  computed  using  heights  at  the  obstacle  horizons. 
As  was  previously  noted  the  values  of  N   for  input  can  range 
from  250  to  400  N-units,  so  it  is  obvious  that  predictions  of 
path  loss  for  anomalous  conditions  are  not  readily  obtained 
using  TR0P0PL0T  which  is  designed  for  1 ong-term  and  medi  an 
input  parameters.   The  effects  of  non-standard  atmospheric 
conditions  along  these  paths  are  considered  in  Section  III. 

After  plotting  the  path  profile  the  required  input 
parameters  were  determined  using  the  methods  described  in 
Appendix  B  and  a  number  of  parameters  were  varied  to  produce 
the  results  described  below. 
2 .    Resul ts 

The  output  of  the  path  profile/routine  for  each 
path  considered  is  shown  in  Figs.  8-19  while  the  corresponding 
input  parameters  used  are  contained  in  Table  I.   All  the  links 
were  analyzed  for  path  loss  at  both  VHF  and  S-band  to 


31 


140 


120 


100 


—    80 

E 


60 


LU 
O0 


O 

CQ 

«c     40 


o 


20     . 


Keyport    (13.11) 


5  10  15 

DISTANCE    (km) 

Figure  8 
Terrain  Profile  for  the  Zelatched  Point  to  Keyport  Path 


32 


300- 


Makah    (365.76) 


200 


Striped    Peak 
(355.4) 


40  60 

DISTANCE    (km) 


Figure   9 
Terrain  Profile   for   the  Makah   to   Striped  Peak  Path 


33 


Lookout    Mtn.     (782.42) 


400 


200 


UJ 

u~> 

> 

o 

CO 

< 
I— 
o 


Keyport    (13.11) 


20  40  60  80 

DISTANCE    (km) 

Figure  10 
Terrain  Profile  for  the  Lookout  Mt.  to  Keyport  Path 


100 


34 


140- 


100- 


LU 
C/7 


o 
en 


o 


4  6 

DISTANCE    (km) 


Figure   11 

Terrain  Profile   for    the   Zelatched  Point   to 
Bangor  Path 


35 


600 


400- 


200 


LU 


o 

CO 


C3 

— 
UJ 


100  150 

DISTANCE    (km) 


Figure   12 

Terrain  Profile  for  the  Mt.  Constitution  to 
Gold  Mt.  Path 


36 


300 


200 


100 


LU 
CO 


o 

ca 


l-H 

LU 


Striped  Peak 
(355.40) 


M t .  Constitution 


40  60 

DISTANCE    (km) 


Figure   13 

Terrain  Profile  for  the  Striped  Peak  to 
Mt.  Constitution  Path 


37 


500 


C.   400 


go     300 


o 

CO 


o     200 


100 


Gold    Mtn.     (536.75) 


Bangor    (  1 2 8-r^t) 2  ) 


5  10  15 

DISTANCE    (km) 

Figure  14 
Terrain  Profile  for  the  Bangor  to  Gold  Mt.  Path 


38 


500 


400 


LU 

<S) 

300 

UJ 

> 

o 

CO 

«=t 

h- 

=c 

CJ3 

1— 1 

200 

100 


Gold  Mtn.  (536.75) 


10         15 

DISTANCE  (km) 


20 


Figure  15 
Terrain  Profile  for  the  Keyport  to  Gold  Mt.  Path 


39 


,ookout    Mtn.    (  .78242) 


OEF  ,    Bangor 
(  .07315) 


20 


40 


60 


80 


100 


DISTANCE    (km) 


Figure   16 
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illustrate  the  increase  in  path  loss  when  frequency  is 
increased.   The  results  are  shown  in  Figs.  20-31.   The  quanti- 
ties (referred  to  in  an  above  sub-section)  which  required 
interpretation  as  to  the  precise  value  to  be  used  in  detailed 
terrain  studies  are  varied  for  a  single  path  so  that  the 
effects  of  this  variation  can  be  shown.   A  sample  path 
(Lookout  Mt.  to  Keyport)  was  selected  for  this  purpose  and 
the  effects  of  several  different  antenna  heights  in  an 
unmodified  version  of  TR0P0PL0T  are  shown  in  Fig.  32.   These 
same  antenna  heights  are  then  entered  into  the  modified 
program  and  the  results  shown  in  Fig.  33.   To  show  the  effects 
of  variations  in  conductivity  and  permittivity  a    and  e  are 
varied  over  the  same  path  with  the  results  shown  in  Fig.  34. 
The  effect  of  variations  in  the  value  of  Ah  is  shown  in 
Fig.  35.   The  quantity  Ah  was  computed  by  taking  the  mean 
value  of  Ah  over  all  the  paths  considered,  Ah,  is  the  tabulated 
estimation  of  terrain  irregularity  found  in  both  Ref.  6  and 
Ref.  7,  Ah   is  the  parameter  computed  by  the  methods  of 
Appendix  B.   This  comparison  is  intended  to  show  the  degree  of 
precision  required  in  choosing  a  parameter  when  uncertainty 
over  the  means  of  selection  is  present. 

The  importance  of  detailed  path  profile  information 
is  pointed  up  in  the  Zelatched  Pt.  to  Keyport  path  (Fig.  8). 
If  TR0P0PL0T  is  used  without  modification,  errors  in  prediction 
can  occur  at  2000  MHz  the  results  of  which  are  shown  in 
Fig.  20.   Note  that  this  output  predicts  free  space  attenuation, 
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when  there  is  clearly  a  large  attenuation  of  the  signal  by 
path  obstacles.   The  same  error  can  occur  when  the  modified 
program  is  used;  primarily  due  to  the  fact  that  the  terrain 
is  represented  statistically.   Large  deviations  from  the 
range  of  terrain  heights  occurring  in  mid-path  may  not  be 
considered  as  an  obstacle  by  the  program  when  in  fact  it  is 
a  major  factor  in  the  attenuation  process.   These  smoothing 
effects  can  be  alleviated  somewhat  by  examining  the  terrain 
profile  and  performing  graphical  analysis,  such  as  is 
described  in  Refs.  1-5,  to  check  the  effectiveness  of  the 
computer  routine  when  its  output  is  questionable.   Another 
possible  source  of  error  can  occur  in  short  paths  where  only 
a  single  knife  edge  obstacle  is  present.   In  this  case  the 
predicted  value  of  path  loss  is  too  high  since  a  non-LOS  path 
loss  is  calculated  for  a  double  knife  edge  case. 
3 .    Conclusions 

It  is  obvious  that  the  quality  of  the  output  of 
this  tropospheric  propagation  prediction  program  depends  on 
which  of  several  possible  assumptions  were  made  concerning 
the  choice  of  input  parameters.   Clearly,  if  information  is 
desired  concerning  a  particular  path  or  path  area  a  path 
profile  should  be  drawn  and  the  ray  path  checked  to  insure 
that  the  values  of  the  corresponding  angles  of  elevation 
(depression)  fall  within  the  program  constraints.   As 
previously  stated,  if  these  values  exceed  the  limits  then 
TR0P0PL0T  must  be  modified  and  a  detailed  terrain  profile 
obtained. 
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The  effect  of  the  variations  of  several  of  the 
parameters  shown  above  allow  some  inferences  to  be  drawn 
concerning  the  care  with,  which  the  input  variables  must  be 
chosen.   It  would  seem  from  the  path  studied  that  variations 
of  a    and  e  have  the  least  effect  on  path  loss  (<1  db),  varia- 
tions in  Ah  have  only  a  small  effect,  particularly  for  longer 
distances,  and  variations  in  antenna  height  have  the  greatest 
effect  on  path  loss.   The  unmodified  version  of  the  program 
appears  more  susceptible  to  these  changes  than  the  version 
modified  to  accept  detailed  profile  information.   The  sensi- 
tivity of  the  transmission  loss  to  antenna  height  value 
variation  was  alluded  to  in  an  earlier  section  and  thus  the 
suggestion  that  care  be  exercised  in  selecting  these  values 
is  well  taken  if  any  degree  of  precision  is  desired.   The 
detailed  computation  of  Ah,  on  the  other  hand,  is  not  necessary. 
By  using  the  tabulated  values  corresponding  to  a  particular 
type  of  terrain  [Ref.  6]  similar  results  can  be  obtained. 

The  transmission  loss  predictions  seem  to  be 
accurate,  or  at  least  within  an  expected  range  of  values  except 
for  short  diffraction  paths.   In  order  to  assess  the  goodness 
of  the  prediction,  comparison  with  empirical  data  for  these 
paths  should  be  accomplished  in  a  later  study. 
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Ill .   PROPAGATION  IN  A  NON-STANDARD  ATMOSPHERE 

The  computer  analysis  conducted  in  the  previous  section 
produced  a  long-term  median  attenuation  value  as  output  and 
consequently  required  that  standard  atmospheric  conditions 
be  assumed.   In  providing  input  under  this  assumption,  the 
"standard"  atmosphere  for  the  locations  under  consideration 
was  characterized  by  the  median  value  of  minimum  monthly 
surface  refractivity  at  sea  level  (N)  corrected  to  emitter 
elevation.   This  element  of  the  prediction  routine  could,  in 
some  cases,  lead  to  an  over-optimistic  expectation  of  link 
performance,  particularly  in  a  locale  where  anomalous  condi- 
tions frequently  occur.   One  method  of  avoiding  this  potential 
source  of  error  is  to  examine  the  statistical  occurrence  of 
these  conditions  and  the  severities  of  the  effect  of  the 
anomalies  on  the  link  in  question.   For  the  purposes  of  this 
study  the  non-standard  conditions  dealt  with  involved  the 
super-refractive,  subrefractive  and  ducting  cases.   The 
statistics  used  in  this  determination  were  derived  from  IREPS 
output  and  the  effects  of  non-standard  refractivity  were 
modeled  by  assuming  certain  values  of  effective  earth's  radius 
for  corresponding  conditions.   Prior  to  the  presentation  of 
these  results  the  theory  involved  in  the  anomalous  propagation 
problem  is  reviewed. 
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A.    THEORY 

The  primary  medium  through  which  electro-magnetic  waves 
with  wavelengths  less  than  1  or  2  meters  propagate  is  the 
troposphere.   This  region,  bounded  below  by  the  earth's 
surface  and  above  by  the  tropopause,  is  approximately  10  km 
thick  and  is  characterized  by  a  general  decrease  in  tempera- 
ture with  height  up  to  a  zone  of  constant  temperature  called 
the  tropopause.   The  tropopause  is  not  a  static  boundary,  but 
has  a  height  which  is  variable  with  both  time  and  latitude. 

The  troposphere  is  usually  assumed  to  be  a  lossless 
dielectric  with  y  =  1  and  a  =  0.   The  index  of  refraction,  n, 
is  then 

n  =  ^r  atm. 

At  the  earth's  surface  n  has  been  found  to  equal  approximately 
1.0003.   Since  air  with  a  higher  water  vapor  content  has  a 
larger  value  of  permittivity,  that  is, 

£wet  air-  £dry  ai  r 

and  since  the  lower  atmosphere  usually  has  the  higher  water 
vapor  content,  the  permittivity  of  the  atmosphere  exhibits  a 
decrease  with  height  to  a  value  of  unity.   This  effect  causes 
the  refractive  index  to  similarly  decrease  with  height.   It 
has  been  ascertained  [Ref.  3]  that  for  temperate  climates  the 
average  variation  near  the  ground  is 


4-r-  =  -  0.039x10"   per  meter 
dh 
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In  order  to  avoid  the  use  of  such  small  numerical  quantities 

a  variable,  N,  designated  as  the  refractivity  or  co- index  of 

refraction,  has  been  defined  as 

6 


N  =  (n-1  )x!0 


N  -  u  n  i  t  s 


and  so  within  one  kilometer  of  the  earth's  surface 


dN 

-7T-  =  -39  N-un  its/ meter 

A  standard  atmosphere  is  then  defined  in  Refs 
fol 1 ows  : 

1  .   £  =  1  ; 
n  =  1 


1 1  and  12  as 


2. 


dT 
dh 


1   °C 


150   m 


3.   N  =  7-^-    (p  +  4810  j) 


p  =  press  ure  (mb ) 

T  =  absol u te  temp .  ( °  K) 

e  =  vapor  pressure  (mb) 


By  assuming  a  linear  refractivity  gradient  it  is  possible 
to  define  an  effective  earth  radius  as  was  done  in  an  earlier 
part  of  this  report.   This  allows  the  rays  representing  the 
radio  wave  to  be  drawn  as  straight  lines  over  a  curved  earth. 
An  alternative  representation  which  enables  the  effect  of 
several  different  values  of  effective  earth  radius  to  be 
shown  on  the  same  profile  is  to  show  the  rays  as  curves  from 
transmitter  to  receiver  over  a  plane  earth.   Details  of 
constructing  this  sort  of  representation  are  presented  in 
Refs .  3  and  1 3. 
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Non-standard  propagation  can  occur  when  the  above  condi- 
tions are  not  met.   This  is  manifested  by  the  manner  in  which 
the  radio  wave  path  is  curved.   Under  normal  conditions  the 
decrease  of  the  refractive  index  with  height  causes  a  downward 
bending  of  the  wave  path.   Changes  in  the  refractive  index 
super-refraction  or  in  certain  cases  an  upward  bending  of  the 
ray  which  is  termed  s ubref racti on .   If  the  refractive  effects 
are  severe  enough,  the  formation  of  ducts  occurs  in  which  the 
majority  of  radio  wave  energy  is  trapped  within  a  narrow 
regi  on  . 

Ducting  can  occur  as  a  result  of  a  number  of  different 
atmospheric  conditions.   The  primary  requirement  in  producing 
this  phenomena  is  that  the  curvature  of  the  radio  wave  must  be 
greater  than  that  of  the  earth.   This  is  produced  by  a  rapid 
change  in  N  with  height  which  is  caused  by  a  sharp  decrease  in 
moisture  with  height  (abnormal  moisture  lapse  rate)  and/or  a 
sharp  increase  in  temperature  with  height  (temperature 
inversion).   The  formation  of  a  ground-  or  surface-based  duct 
can  occur  when  warm  dry  air  flows  from  land  over  water. 
Moisture  evaporates  from  the  water  into  the  lower  layers  of 
air  cooling  the  air  and  producing  an  increased  moisture  lapse 
rate  and  a  temperature  inversion.   The  downward  motion  of  warm 
dry  air  is  frequently  associated  with  the  clear  weather  found 
on  the  eastern  side  of  vast  high  pressure  regions  in  the  lower 
and  middle  atmosphere.   Subsidence  such  as  this  can  also 
result  in  elevated  ducts  which  are  found  to  be  strongest  and 
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lowest  in  fair  weather  and  highest  and  weakest  near  storms. 
Surface  ducts  are  also  formed  in  certain  regions  by  nocturnal 
cooling  but  unlike  the  evaporation  duct  this  sort  of  anomaly 
if  found  over  land,  particularly  desert  regions. 

The  opposite  effect  to  super-refractive  ducting  is  the 
upward  bending  of  radio  waves  caused  by  an  increase  in 
refractivity  with  height.   This  condition  typically  occurs 
in  moist  and  cloudy  regions. 

The  determining  factor  in  all  the  above  cases  is  the 
refractivity  gradient.   Tabl e  II  bel ow ,  found  in  Ref.  11,  shows 
the  range  of  refractivity  gradients  over  which  anomalous 
conditions  prevail. 


Table  II.   Range  of  refractive  index  gradient  for  differing 
types  of  propagation  conditions. 


p 

(mb) 

h 
(km) 

Sub- 
Refractive 

Unstratified 

Super- 
Refract  i  ve 

1000-850 

0-1 .46 

-dn/dh^O 

20£-dn/dh^60 

100<_-dn/dh 

850-700 

1 .46-3.01 

-dn/dh^O 

20£-dn/dh£50 

80^-dn/dh 

700-600 

3.01 -4.20 

-dh/dh<_0 

20£-dn/dh^40 

70^-dn/dh 

600-500 

4.20-5.57 

-dh/dh<0 

20£-dn/dh<_30 

50<_-dn/dh 

500-400 

5.57-7.  18 

-dn/dh^0 

20£-dn/dh^25 

40^-dn/dh 

The  refractivity  gradient  can  be  related  to  the 
effective  earth  radius  and  the  resulting  ray  trace  when 
drawn  on  the  terrain  profile  corresponding  to  this  radius 
will  show  the  effect  of  the  anomaly  on  a  particular  path. 
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The  following  holds  for  most  conditions  [Ref.  3] 


effective  radius  =  k  = 


1 


1+6.4x10   dN/dh 


From  this  relationship  it  can  be  seen  that  for  the  standard 
refractivity  gradient,  which  was  defined  as  -39  N-un its/km, 
the  effective  earth  radius  would  be  equal  to  4/3.   The 
presence  of  anomalous  conditions  thus  cause  an  effective 
earth  radius  greater  or  less  than  4/3.   The  interrelationship 
of  the  refractivity  gradient,  the  radius  of  curvature  of  the 
wavefront,  the  effective  radius,  and  the  assoicated  types  of 
refraction  are  shown  in  Table  III  [Ref.  3]. 

Note  that  there  is  some  difference  in  the  exact  gradient 
value  which  defines  s uper- re f racti on  .   The  s uper-ref racti vi ty 
gradient  referred  to  in  Table  II  includes  both  the  extended 
range  and  the  ducting  conditions,  while  that  referred  to  in 
Table  III  includes  only  the  ducting  case. 

In  the  actual  paths  under  study  several  values  of 
effective  earth  radius  were  used  to  simulate  the  effect  of 
anomalous  conditions.   Table  IV  [Ref.  15]  provides  a  guide 
to  the  use  of  k  as  an  estimate  of  propagation  conditions  for 
99 . 9%-99 .99%  path  reliability. 
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Table  IV.   K  Factor  Guide 


Propagation  Conditions 


Perfect 


Ideal 


Average 


Difficult 


Bad 


Weather 


Standard 
atmosphere 


No  surface 
layer  or 
fog 


Substandard 
light  fog 


Surface 
1 ayers  , 
ground 
fog 


Fog 

moi  s  ture 

over 

water 


Typi  cal 


Temperate 
zone  ,  no 
fog  ,  no 
ducti  ng  , 
good  atmos 
p  h  e  r  i  c  mix 
day  and 
ni  ght 


Dry,  moun 
tai  nous  , 
no  fog 


Flat, 

temperate , 
some  fog 


Coas  tal 


Coas  tal  , 
water , 
tropical 


K  Factor 


1  .33 


1-1  .33 


0.66-1  .0 


0.66-0.5 


0.5-0.4 


The  path  profiles  from  Keyport  to  Gold  Mt.,  Keyport  to 
Lookout  Mt.,  and  Bangor  to  Lookout  Mt.  are  shown  for  two 
values  of  K;  0.5  which  represents  the  difficult  case  as  either 
when  one  antenna  is  inside  a  duct  and  the  other  outside,  or 
the  subref racti ve  condition,  and  -3.57  which  represents  a 
sub  refractive  gradient  of  -200  N-un its/km.   The  effect  of  the 
change  in  earth  radius  shown  in  Figs.  36-41  is  most  noticeable 
in  the  Keyport  to  Lookout  Mt.  path.   Using  graphical  methods 
of  analysis  found  in  Refs.  1  and  3  the  approximate  degradation 
in  link  performance  can  be  estimated.   While  little  change  is 
observed  in  the  Keyport  to  Gold  Mt .  link  due  to  the  shortness 
of  the  path  length,  the  Keyport  to  Lookout  Mt.  path  with  K  =  h 
exhibits  an  additional  loss  of  up  to  60  db  over  that  found  in 
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Figure   36 
Path  Profile   for  Lookout  Mt.    to  Keyport    (K=l/2) 
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Figure  37 
Path  Profile  for  Lookout  Mt.  to  Keyport  (K=l) 
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Figure  38 
Path  Profile  for  Keyport  to  Gold  Mt.  (K=l/2) 
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Figure  39 
Path  Profile  for  Keyport  to  Gold  Mt.  (K=l) 
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Figure  40 
Path  Profile  for  Lookout  Mt.  to  Bangor  (K=l/2) 
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Path  Profile  Lookout  Mt.    to   Bangor    (K=l) 


100 


76 


standard  atmospheric  conditions.   This  same  path  however 
suffers  no  apparent  degradation  for  K  =  1  as  the  Fresnel  zone 
clearance  remains  greater  than   .6.   The  Lookout  Mt.  to 
Bangor  path  degradation  was  found  to  be  12-30  db  for  K  =  h. 
and  6-15  db  for  K  =  1.   From  these  results  it  can  be  seen 
that  changes  in  effective  earth  radius  can  have  a  significant 
effect  on  path  loss,  particularly  as  the  path  length  increases 

The  profile  representation  using  N  data  is  the  most 
commonly  used  in  the  literature,  however  the  depiction  of 
ducting  or  trapping  types  of  profiles  is  often  done  with 
modified  or  M-unit  profiles  where 


M  =  N  +  i|l  x  10 
h 


a  =  actual  earth  radius 
h  =  height  above  sea  level 


This  adjustment  to  the  N-profile  is  used  in  the  anomalous  case 
since  straight  rays  above  a  curved  earth  now  become  curved 
rays  above  a  planar  earth  and  the  duct  can  be  more  clearly 
shown.   The  variability  of  M  curves  corresponding  to  a 
particular  propagation  condition  is  shown  in  Table  V  which 
was  derived  from  Ref.  16. 

The  M-curve  is  thus  a  transformation  in  which  the 
relative  curvature  between  the  normal  of  the  electro-magnetic 
wavefront  and  the  surface  of  the  earth  is  unchanged.   A  duct 
is  formed  whenever  the  M  curve  has  a  relative  minimum.   The 
base  of  an  elevated  duct  is  the  height  on  the  M  curve  below 
the  minimum  that  has  the  same  M  value  as  the  relative  minimum. 
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When  no  such  height  exists,  the  base  of  the  duct  is  at  the 
ground  and  a  so-called  ground  based  duct  exists  [Ref.  12]. 
The  main  disadvantage  in  the  use  of  M  data  is  that  above  the 
duct  N  is  grossly  overcorrected . 


Table  V.   Variability  of  M  curves  with  type  of  transmission 
conditions . 


Increases  linearly 


M-Curve 

Straight  line,  with 
positive  slope 

Slope  decreases  near 
earth  surface  with 
rays  curving  up 

Slope  increases  near 
surface  with  less 
upward  curve 

Vertical,  no  curve 


Type  of 
Transmission 

Standard 


Subs  tandard 


Supers tandard 


Very    great 
coverage 


Decreases  wi  th 
height 


SI  ope  linear  and 
negative  rays 
curve  down 


Trapping  or 
ducting 


Having  demonstrated  the  effect  of  anomalous  conditions 
on  several  propagation  paths  and  considered  the  theory 
involved  the  next  step  is  to  determine  the  frequency  with 
which  these  conditions  occur.   The  method  chosen  to  accomplish 
this  was  to  access  the  IREPS  historical  file  for  the  geograph- 
ical area  in  question  as  described  in  the  following  section. 


B.    INTEGRATED  REFRACTIVE  EFFECTS  PREDICTION  SYSTEM  (IREPST 

IREPS  is  an  experimental  system  ultimately  designed  to 
provide  the  capability  for  on-board  assessment  of  the  effect 
of  atmospheric  anomalies  on  sensor  performance.   The  hardware 
portion  of  the  system  consists  of  a  mi ni -computer  with  an 
interactive  graphics  terminal  on  which  the  required  inputs  are 
entered  and  the  results  displayed  in  an  easily  understood 
graphic  format.   The  information  stored  in  the  systems  memory 
includes  a  refractivity  library  containing  long-term  meteoro- 
logical statistics  on  the  occurrence  of  ducting  conditions  as 
a  function  of  location.   This  library  data  is  obtained  from 
radiosonde  stations  throughout  the  world  and  can  be  augmented 
by  on-scene  re f ractome ter  or  radiosonde  data  input  in  order 
to  tailor  the  output  to  a  particular  location.   In  the  case 
of  the  Washington  paths,  there  were  only  two  upper  air  sounding 
stations  in  the  vicinity  of  the  sites.   One  is  located  at 
Tatoosh  Island  (48°-24'N,  124°-42'W)  and  the  other  at 
Quillayute  (48°-00'N,  124°-36'W).   The  historical  propagation 
summaries  for  each  location  are  shown  in  Figs.  42  and  43. 
Note  that  while  the  statistics  generated  are  geared  toward 
radar  performance,  the  frequencies  to  be  used  in  the  link  are 
included  so  that  some  information  concerning  potential 
anomalies  can  be  derived.   The  historical  meteorological  data 
for  each  station,  shown  in  Figs.  44  and  45,  provides  informa- 
tion concerning  median  conditions  in  the  vicinity  of  the 


2  F  o  r  a  detailed  description  of  the  system  see  Ref.  10 
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sounding  station  and  shows  the  seasonal  and  diurnal  variation 
in  this  data.   Based  on  this  data,  the  program  then  can 
generate  a  series  of  profiles.   Figures  46  and  47  are  M  and 
N  profiles  in  the  vicinity  of  the  two  sounding  stations. 
Also  available  in  the  program  output  capabilities  is  a  ray 
trace  diagram  from  a  particular  emitter.   This  is  useful  in 
a  qualitative  sense  in  that  it  shows  potential  areas  of  non- 
coverage  which  may  be  important  in  maintaining  a  ground  to  air 
communications  link.   The  first  set  of  traces  demonstrates  the 
effect  on  the  ray  trace  diagram  of  selected  sites  for  the 
median  ducts  found  in  the  Tatoosh  Island  profile  and  are 
shown  in  Figs.  48-52.   A  second  set  of  traces  are  shown  in 
Figs.  54-57  exhibiting  the  effect  of  an  elevated  duct  between 
1100  and  1400  feet.   The  profile  responsible  for  this  condi- 
tion is  shown  in  Fig.  53.   The  third  set  of  traces  show  a 
duct  between  2300  and  2900  feet  in  Figs.  59  and  60  and  the 
corresponding  profile  in  Fig.  58. 

Besides  the  above  mentioned  output  format  the  IREPS 
program  can  also  provide  path  loss  information,  with  a 
detection  threshold  for  a  particular  sensor  superimposed  on 
the  distance  vs.  dB  display  and  also  a  coverage  diagram  for 
a  particular  sensor.   These  two  output  diagrams  were  not 
produced  for  this  study  but  an  example  of  each  is  included  in 
Figs.  61  and  62  [Ref.  10]. 
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This  system  is  not  yet  operational  and  is  still  under- 
going testing.   The  value  of  the  output  formats  in  evaluating 
link  anomalous  propagation  potential  is  clearly  obvious.   If 
data  were  available  for  each  of  the  sites  analyzed  in  the 
study  a  more  exact  picture  of  the  statistics  of  propagation 
would  be  available. 

Based  on  the  available  data  for  the  Washington  area 
several  inferences  can  be  made  concerning  anomalous  propaga- 
tion.  For  the  frequencies  in  question  the  probability  of 
extended  ranges  or  ducting  is  between  8  and  10  percent  while 
no  statistics  are  presented  for  be! ow  normal  refraction.   For 
the  locations  from  which  the  soundings  are  available  it  can 
be  seen  that  surface  ducting  occurs  most  often  in  the  daytime 
during  Spring  and  that  the  ducts  are  predominantly  less  than 
12  meters  in  thickness.   The  lowest  occurrence  of  surface 
ducts  at  Quillayute  is  in  the  spring  night  which  also  has  the 
highest  incidence  of  elevated  ducts.   At  Tatoosh  which  is  a 
more  classically  marine  climate,  there  are  no  ducts  observed 
to  occur  during  the  autumn  and  winter  either  at  night  or 
during  the  day.   The  differences  in  the  median  refractivity 
at  sea  level  are  not  significant  and  both  values,  332  for 
Tatoosh  and  324  for  Quillayute,  would  produce  an  effective 
earth  radius  greater  than  4/3,  and  consequently  better  than 
average  propagation  conditions  prevail.   Until  meteorological 
data  is  available  for  the  actual  sites  the  ray  trace  diagrams 
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can  provide  only  an  indication  of  the  effect  of  a  certain 
profile  on  the  emitter  in  question  rather  than  actual  informa- 
tion concerning  real  conditions  present  in  the  link  components 


103 


APPENDIX  A 

Path  Plotting  Program 

In  order  to  plot  a  path  profile  on  a  curved  earth  surface 
either  commercially  available  profile  paper  must  be  used  or 
the  profile  can  be  manually  or  computer  generated.   Instruc- 
tions for  the  construction  of  path  profile  paper  for  various 
values  of  earth  radius  can  be  found  in  Ref.  13.   Computer 
generated  profiles  can  be  easily  produced  if  a  plotting 
package  is  available.   For  the  purposes  of  this  study  a 
subroutine  known  as  DRAWP  was  locally  available.   Other 
alternative  plotting  methods  included  the  TEKTRONIX  4012  and 
HP  9830  plotting  packages,  both  of  which  are  available  at  the 
Naval  Postgraduate  School  . 

The  program  used  made  use  of  a  transformation  which 
allowed  the  terrain  profile  to  be  plotted  on  a  linear  graph 
by  adjusting  the  height  at  a  particular  distance  for  the 
effective  curvature  of  the  earth.   The  profile  is  drawn  by 
the  CALCOMP  plotter  along  with  the  earth  surface  at  sea  level. 
The  program  also  fits  a  line  to  the  terrain  data  by  the  method 
of  least  squares  and  forms  the  difference  between  the  value 
of  height  found  on  the  line  and  the  actual  value  of  the 
terrain  height  at  the  same  distance.   Aside  from  the  actual 
terrain  data  points  the  only  input  parameters  necessary  to 
obtain  this  output  are  either  the  refractivity  gradient  or 
the  surface  refractivity,  and  the  number  of  points  input. 
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The  output  provided  by  the  program  includes  a  plot  of  the 
terrain  on  the  effective  earth  surface,  a  table  of  height  vs 
input  distance,  a  list  of  the  terrain  heights  in  kilometers, 
the  effective  earth  radius,  and  the  difference  between  each 
point  on  a  least  square  curve  and  the  actual  terrain  data 
points  without  the  correction  for  earth  curvature  applied. 
In  order  to  perform  the  curve  to  fit  terrain  referred  to  in 
Appendix  B,  a  new  set  of  points  must  be  used  taking  only 
those  the  points  visible  to  both  transmitter  and  receiver  instead  of 
the  entire  terrain  profile. 

The  format  of  the  input  data  deck  is  presented  below: 
First  Data  Card: 

Column  1:   Code  "1".  surface  refractivity  will  be 

'used  as  input 

Code  "2".  refracti  vi  ty  gradient  will  be 
used  as  input 


Second  Data  Card: 

Column  1-10:   Number  of  data  points 

(right  justified  integer  <_  200) 

Column  10-20:  Refractivity  or  gradient 
(right  justified  integer) 


Third  Through  Last  Data  Card: 

The  distance  points  in  kilometers  are  listed  in  order 
usinq  real  numbers  (i.e.,  F10.5  format)  in  columns 
1-10   11-20,  ...,  61-70.   The  heights  above  sea  level 
in  feet  follow  directly  after  the  distance  listing 
in  the  same  format.   A  new  card  is  not  started  unless 
the  last  distance  point  was  listed  in  Columns  61-70 

An  example  of  the  latter  cards  is  shown  below 
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For  the  table  of  distance  vs.  height  bel 


ow 


Distance 
(km) 

Height 
(feet) 

0.0 

420. 

1.0 

500. 

2.0 

10. 

3.0 

125.5 

4.0 

780. 

5.0 

1000. 

The  final   data  cards  would  be  as  follows: 
Column:        1-10       11-20       21-30       31-40       41-50 


51-60       61-70 


Card  3  0.0         1.0  2.0  3.0  4.0 

Card  4       500.  10.         125.5       780.       1000. 


5.0 


420. 
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APPENDIX  B 
INPUT  PARAMETER  COMPUTATION 

A.  DETERMINATION  OF  DELTA-H 

The  program  described  in  Appendix  A  was  used  to  produce 
from  the  path  geometry  the  parameters  d.  ,  ~>  9  i  ?  anc*  An* 
The  computation  of  Ah  is  only  partially  computer-aided  by  the 
fitting  of  a  straight  line  to  the  uncorrected  profile  heights 
and  subtracting  the  height  of  the  line  from  the  corresponding 
terrain  height.   This  difference  is  output  in  the  array  D(I). 
The  interdecile  range  is  then  extracted  manually  from  this 
data  by  finding  the  value  of  the  lower  boundary  of  the  90th 
percentile  and  subtracting  from  it  the  value  of  the  upper 
boundary  of  the  10th  percentile.   An  example  of  the  fit  of 
this  curve  to  actual  terrain  data  is  shown  in  the  main  body  of 
this  study. 

B.  DETERMINATION  OF  HE!  ,2 

For  line  of  sight  paths  the  following  determination  was 
made  for  effective  antenna  heights: 

1.  A  straight  line  is  fitted  to  the  corrected 
terrain  profile  heights. 

2.  The  line  is  defined  by 

B(I)  =  AO  +  A1(X(I)) 

in  the  program  defined  in  Appendix  A  and  should 
be  modified  for  earth  curvature  by 

B'(D  =  B(I)  -  X2/2A 
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3.  The  points  chosen  in  the  actual  profile  to  which 
this  line  is  fitted  should  be  only  those  points 
which  are  visible  to  both  transmitter  and  receiver. 

4.  The  effective  antenna  height  is  then  the  difference 
between  the  fitted  line  and  the  height  of  the 
antenna  above  sea  level. 

The  procedure  described  above  defines  a  reflecting  plane 

between  the  transmitter  and  receiver. 

If,  on  inspection  of  the  actual  terrain  and  the  least 
square  curve,  a  good  fit  does  not  appear  to  have  been  obtained 
due  to  the  roughness  of  the  terrain  other  methods  of  obtaining 
effective  antenna  heights  are  available.   If  the  foreground 
of  the  antenna  represents  a  good  reflecting  surface  then  the 
antenna  height  over  ground  can  be  used,  or  multiple  least 
square  curves  can  be  used  if  the  path  consists  of  a  number  of 
reflecting  planes.   In  the  event  that  terrain  characteristics 
defy  this  sort  of  analysis,  curves  and  formulae  are  available 
in  Annex  III  to  Ref.  5  which  provide  alternate  methods  for 
computing  effective  height. 

The  curve  to  fit  terrain  is  also  used  in  the  determination 
of  effective  heights  for  the  case  of  knife  edge  diffraction 
which  is  assumed  to  occur  when  the  path  length  is  less  than 
the  smooth-earth  radio  horizon  for  each  antenna.   In  this  case 
two  curves  are  fitted  to  each  of  the  diffracted  paths  forming 
a  reflecting  plane  for  each  antenna. 

For  trans-horizon  paths  the  effective  heights  are  obtained 
by  considering  the  actual  height  of  the  antenna  above  sea  level 
and  subtracting  from  that  height  the  average  height  of  the 
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terrain  between  the  transmitter  and  receiver.   To  illustrate 
this,  consider  a  link  with  transmitter  height  above  sea  level 
hto'   The  nei9hts  at  N  equi-distant  points  are  selected  and 
the  mean  of  the  central  80%  of  these  values  is  computed  to 


produce  h.  where 


.9N 

h"t  =  TW  £  hti 
i  =  .  1  N 


e.g.,   N  =  3 1  ;  h 


27 

2T   £  hti 
i  =  3 


i  =  0,1  ,2, 


,30 


The  effective  transmitter  height  is  then 


hfll  =  h.   -  h.     for   h.  <  h. 
el     to    t  t  —  to 


If  the  mean  value  is  greater  than  the  height  above  sea 
level  then  the  structural  height  of  the  antenna  above  ground 
is  used. 

C.    DETERMINATION  OF  DL1,2  and  TE1,2 

The  angular  distance  9  is  readily  obtained  from  the 
geometry  of  the  profile  as  shown  in  Fig.  2.   9  ,  ,  and  8  « 
are  measured  and  8   is  calculated  in  the  program  as  the 
maximum  of  either  9  ,  +  9  ~    or  -d^/a.   The  distance  dL  is  the 
sum  of  the  distances  to  the  obstacle  horizons,  d.^  and  dL2 
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D.    MODIFICATION  TO  TROPOPLOT 

With  the  above  means  of  determining  the  parameters 

9  -,    o,  d,  ,  o,  TROPOPLOT  was  modified  as  follows: 
el  ,2    11,2' 

1.  Delete  the  following  statements  from  the  main 
program : 

IF  (H1G.LE.2.))  GO  TO  12 

IF  (H1G.GE.2.) .0R.H1G. LE.5. ) )  GO  TO  10 

Zl  =  5.0 
GO  TO  11 

10  Zl  =  1 .0+DSIN(3.1415927*HlG/10.0) 

11  HIE  =  H1G+Z1*DEXP(-2.0*H1G/DH) 
IF  (H2G.LE.2.0)  GO  TO  15 

IF  (H2G.GE.2.0.0R.H2G.LE.5.0)  GO  TO  13 
Z2  =  5.0 
GO  TO  14 

13  Z2  =  1  .0  +  DSIN(3.1415927*H2G/10.0) 

14  H2E  =  H2G+Z2*DEXP(-2.0*H2G/DH) 

15  CONTINUE 

DL1  =  DLS1*DEXP(-.07*DSQRT(DH/DMAX1(5.0,H1E))) 
DL2  =  DLS2*DEXP(-.07*DSQRT(DH/DMAX1(5.0,H2E))) 
TE1  =  (  .00065/DLS1  )*(  (  DLS1 /DL1 -1  .  )  *  DH  -  3  .  0  7  7*H  1  E) 
TE2  =  (.00065/DLS2)*((DLS2/DL2-l .  )  *DH -3 . 077*H2E ) 

2.  Insert  the  following  statements 

909  FORMAT  (4F10.5) 

READ(5,909)  DL1 ,  DL2,  TE1  ,  TE2 
after  the  statement 
22  FORMAT  (4F10.5) . 

Having  accomplished  these  insertions  and  deletions  the  input 
card  deck  as  shown  in  Ref.  7  is  changed  as  shown  below  for 
user-supplied  data  only: 
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First  Data  Card: 

Column  01-10   "Type  of  Terrain"  -  Integer 

Code:  4   User  Terrain  Data  Used 
No  data  suppression 

5   User  Terrain  Data  Used 

No  "output  parameters"  printed 

Column  11-20   "DB  Loss"  -  Integer 

Code:  0   No  DB  data  desired 

(omit  2nd  data  card) 

1   DB  data  desired 

Column  21-30   "Distance  Between  Antennas"  -  Real 
(1-2000)  KM 


Second  Data  Card: 

Column  01-10   "TX  Power  Out"  (Watts) 
Positive  Real  Number 


Column  11-20 

Column  21-30 

Column  31-40 

Column  41-50 


"TX  Antenna  Gain"  (db) 
Positive  Real  Number 

"RX  Antenna  Gain"  (db) 
Positive  Real  Number 

"Transmitter  Line  Loss"  (db) 
Positive  Real  Number 

"Required  Receiver  Line  Loss"  (db) 
Positive  Real  Number 


Column  51-60   "Receiver  Sensitivity"  (dbm) 
Positive  Real  Number 
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Third    Data    Card: 

Column    01-10      "Surface    Ref racti vi ty " 

Integer    (250-400)    N-units 

Column  11-20   "Surface  Conductivity" 
Real  (mhos/meter) 

Column  21-30   "Relative  Dielectric  Constant" 
Real 

Column  31-40   "Interdecile  Range" 
Del ta-H 


Fourth  Data  Card: 

Column  01-10   "Antenna  Polarization" 
Real 

Code:  01.00  Vertical  Polarization 
Code:  -1.00  Horizontal  Polarization 


Column  11-20 
Column  11-20 

Column  31-40 


"Frequency"  Real 
(20-40,000) 

"Calculated  Antenna  Height  of 
Transmitting  Antenna" 
Real  (0.5-3000.0)  Meters 

"Calculated  Antenna  Height 
of  Receiving  Antenna" 
Real  (0.5-3000.0)  Meters 


Fifth  Data  Card: 

Column  01-10   "Transmitter  Obstacle  Horizon,  DL1  "  (km) 
Real 

Column  11-20   "Receiver  Obstacle  Horizon,  DL2"  (km) 
Real 

Column  21-30   "Transmitter  Elevation  (Depression) 
Angle"  (radians) 
Real  Positive  (Negative) 

Column  31-40   "Receiver  Elevation  (Depression) 
Angle  TE2"  (radians) 
Real  Positive  (Negative) 
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